Neutron skin of 208 Pb, nuclear symmetry energy, and the parity radius experiment 
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A precise determination of the neutron skin Ar np of a heavy nucleus sets a basic constraint on 
the nuclear symmetry energy (Ar np is the difference of the neutron and proton rms radii of the 
nucleus). The parity radius experiment (PREX) may achieve it by electroweak parity- violating 
electron scattering (PVES) on 208 Pb. We investigate PVES in nuclear mean field approach to 
allow the accurate extraction of Ar np of 208 Pb from the parity- violating asymmetry A pv probed 
in the experiment. We demonstrate a high linear correlation between A pv and Ar np in successful 
mean field forces as the best means to constrain the neutron skin of 208 Pb from PREX, without 
assumptions on the neutron density shape. Continuation of the experiment with higher precision 
in A pv is motivated since the present method can support it to constrain the density slope of the 
nuclear symmetry energy to new accuracy. 
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New interest in masses and density distributions of nu- 
clei is being prompted by the production of rare isotopes 
in radioactive beam facilities [1 . Exciting phenomena 
discovered in these isotopes such as thick skins, halos, 
and new shell closures urge better understanding of neu- 
trons in nuclei. Yet, our knowledge of neutron density 
distributions is limited even in the stable nuclei. As neu- 
trons are uncharged, neutron densities have been probed 
mostly by nucleon scattering [2 , 3j, a scattering [4 , and 
nuclear effects in exotic atoms [5j [6]. Even if some of 
these experiments reach small errors, all hadronic probes 
require model assumptions to deal with the strong force 
introducing possible systematic uncertainties. 

Parity- violating electron scattering (PVES) was sug- 
gested as a model-independent probe of neutron densities 
[7]. An electroweak probe is not hindered by the com- 
plexity of the strong force and the reaction mechanism 
with the nucleus needs not be modeled [THD], similarly 
to clean electron scattering for nuclear charge densities. 
The novel parity radius experiment (PREX) at the Jef- 
ferson Lab [9[ [10] aims to measure the parity- violating 
asymmetry A pv in polarized electron scattering on 208 Pb 
to 3% accuracy. This accuracy is estimated to constrain 
the neutron rms radius r n of 208 Pb to 1% [10]. Cur- 
rently, r n of 208 Pb is uncertain by ~2% and data may be 
model dependent [2-6, 11 ; in contrast, the charge radius 
of 208 Pb is accurately known as r ch = 5.5010(9) fm [12]. 
In recent years it has been established that the neutron 
skin thickness Ar np = r n — r p (difference of the neutron 
and proton rms radii) of 208 Pb is strongly correlated with 
the density dependence of the nuclear symmetry energy 
around saturation [T3HT7] . Knowledge of the density de- 
pendence of the nuclear symmetry energy is a cornerstone 
for drip lines, masses, densities, and collective excitations 
of neutron-rich nuclei [T3HT9] , flows and multifragmenta- 
tion in heavy- ion collisions [201 El] , and for astrophysi- 
cal phenomena like supernovae, neutrino emission, and 



neutron stars p~5] [22H24] . A constraint from PREX on 
Ar np of 208 Pb is thus regarded as a landmark for isospin 
physics. In addition to being important for its own sake, 
it has broad implications for different communities of nu- 
clear physics and astrophysics. Fostered by the seminal 
study of Ref. [9 , PREX completed an initial run in 2010. 
First analyses [10] show the validity of the experimental 
technique, the adequacy of instruments, and that sys- 
tematic errors are under control. Additional beam time 
is now under request to attain the planned 3% accuracy 
in the parity- violating asymmetry A pv [ID] . 

The direct output of PREX is the value of the asymme- 
try A pv at a single scattering angle [9j [ID] . The neutron 
rms radius r n of the nucleus may be deduced only if a 
shape for the neutron density such as a two-parameter 
Fermi function [9J is assumed. A systematic uncertainty 
in the analysis is unavoidable in this way. Here, we pro- 
vide a different and accurate strategy to deduce r n and 
Ar np from PREX that removes this problem. By study 
of PVES on 208 Pb in successful nuclear mean field (MF) 
forces of wide use in nuclear research and astrophysical 
applications, we reveal a high linear relation between 
Ar np and A pv that allows one to extract r n and Ar np 
from A pv model and shape independently. Moreover, 
our approach unifies the extraction of Ar np from A pv 
with the same framework where Ar np is correlated to the 
symmetry energy. We show that the present method can 
support PREX to narrow down the value of the density 
slope of the nuclear symmetry energy to novel accuracy. 
This result provides a new and important motivation to 
continue the experiment to increased precision. 

Electrons interact with nuclei by exchanging photons 
and Z° bosons. The former mainly couple to pro- 
tons and the latter to neutrons because, opposite to 
the nucleon electric charges, the neutron weak charge 
Qyy = — 1 is much larger than the proton weak charge 
Qw = 1 - 4 sin 2 W « 0.075 (# w being the Weinberg an- 
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gle). Therefore, electron scattering can probe both the 
electric and the weak charge distributions in a nucleus [7]- 
[9]. PREX measures the elastic differential cross sections 
da±/d£t for incident electrons of positive or negative he- 
licity. The parity- violating asymmetry, 

pv ~\dn cm)/ \dn + dn) [ ' 

for massless electrons (it is m e /p e ~ 0.0005 at PREX 
energy), is sensitive to the parity- violating term induced 
by the weak interaction in the scattering amplitude. Ac- 
cording to their helicity, electrons interact with a po- 
tential ^coulomb (r) ± GpPw(j)/2 3 / 2 , with Gp the Fermi 
constant and pw the weak density of the target [7H9]. 
We solve the associated Dirac equation via the exact 
phase-shift analysis in distorted-wave Born approxima- 
tion (DWBA) [11] to compute A pv . Our benchmarks 
are the pointlike densities of protons p p (r) and neu- 
trons p n (r) calculated self-consistently in MF models. 
We fold p p (r) and p n (r) with electromagnetic proton 
and neutron form factors to obtain the charge density 
[IT], and with electric form factors for the coupling to 
a Z° to obtain the weak density [9j [TTJ [25] : pw(r) — 
Jdr'{4G?(r') N Pn (\r - r'|) + 4Gf (r') Z Pp (\r - r'|) }. 

Though not useful for realistic calculations, it is worth 
recalling the Born approximation (B A) to A pv pj [9] : 

as it nicely illustrates that A pv relates to the neutron 
and proton nuclear form factors F n ^ p (q). Furnstahl [T4] 
showed that F n (q) = (47r) _1 jd 3 r jo(qr)p n (r) is at low 
momentum transfer q strongly correlated with r n of 
208 Pb in nuclear MF models, evidencing that PREX 
would directly constrain the neutron radius and the sym- 
metry energy. Realistic DWBA calculations of A pv in MF 
models can be found in EJ [25] [26] . 

At the optimal kinematics of PREX the electron beam 
energy is 1.06 GeV and the scattering angle is 5° (#i a b~ 
0.47 fin -1 ) [10 . We compute A pv in DWBA at this 
kinematics in a comprehensive large sample of 47 nu- 
clear MF interactions. We display the results in Fig. [I] 
as a function of the neutron rms radius of 208 Pb. To 
prevent eventual biases in our study, we avoid including 
more than two models of the same kind fitted by the 
same authors and protocol. We also avoid models yield- 
ing a charge radius of 208 Pb away from experiment [12] by 
more than 1% (same level as the 1% pursued by PREX in 
r n ). The considered models rest on very different theo- 
retical grounds, from nonrelativistic models of zero range 
(models HFB, v090, and those starting with S or M) 
or finite range (D1S, DIN, BCP), to relativistic models 
with meson self-interactions (NL and PK models, FSUG- 
old, Gl, G2, TM1), density-dependent vertices (DD-ME, 
RHF-PK) or point couplings (DD-PC1, PC-PK1, PC- 
PF1) [IIJ[I7HI9]. (NL3.s25 and PKl.s24 are variants of 
NL3 and PK1 giving Ar np = 0.25 and 0.24 fm in 208 Pb.) 
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FIG. 1: Parity- violating asymmetry for 208 Pb at the kine- 
matics of PREX against the neutron radius of 208 Pb in nu- 
clear models. The linear fit is 10 7 A pv = 25.83 — 3.31r n - The 
inner/outer colored regions depict the loci of the 95% con- 
fidence/prediction bands of the regression (see e.g. Ch. 3 of 
[27]). An assumed sample measurement A pv = 0.715 ppm of 
3% accuracy and its projection on the r n axis are also drawn. 

All such models accurately describe general properties of 
nuclei such as binding energies and charge radii along the 
periodic table. However, one readily sees in Fig. [I] that 
the predicted r n of 208 Pb varies largely, from 5.55 to 5.8 
fm, as the isovector channel of the nuclear models is little 
constrained by current phenomenology. The models with 
softer (stiffer) symmetry energy at saturation density [TT] 
yield smaller (larger) r n and larger (smaller) A pv . One 
notes that the information encoded in the models implies 
a value of about 0.67 to 0.75 ppm for A pv at PREX kine- 
matics. A significant linear trend is found between A pv 
and r n (the correlation coefficient is r = 0.974). 

As the experimental value of A pv is not yet avail- 
able, we have chosen for study a plausible test value 
A pv = 0.715 ppm of 3% accuracy, depicted in Fig. HI The 
assumed sample measurement of A pv determines through 
the linear fit shown in Fig. [T] a fiducial neutron rms ra- 
dius r n = 5.644 ±0.065 fm, within typical values deduced 
from hadronic probes [2-6 . Note that a 3% accuracy in 
A pv does lead to ~1% accuracy in r n , thereby supporting 
the expectations of PREX. It is to be pointed out that 
the analysis described in this paper is actually indepen- 
dent of the exact value of the parity- violating asymmetry. 
Thus, once the experimental value of PREX is known, 
one can repeat the same type of analysis using the actual 
A pv instead of our test value. We also plot in Fig. [T] the 
confidence band of the regression (boundary of the possi- 
ble straight lines) and the so-called prediction band (the 
wider band that basically coincides with the envelope of 
the models in the figure) at 95% confidence level [27] . 

While one first thinks of using a PREX extraction of 
r n to constrain Ar np of 208 Pb, we show in Fig. that 
A pv and Ar np have themselves a very high linear depen- 
dence (the correlation coefficient is 0.995). The small 
fluctuation of A pv with the charge density is more ef- 
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FIG. 2: Same as Fig. [l] against the neutron skin of 208 Pb. 
The linear fit is 10 7 A pv = 7.88 — 3.75Ar np . The correlation is 
found to be quite stable: for example, if we remove the forces 
excluded by the depicted test constraint, then r = 0.990. 
The figure also shows the points calculated with the neutron 
densities deduced from experiment in Refs. [21 El E] - 



fectively removed by analyzing A pv vs r n — r p . Ac- 
tually, the correlation of A pv and Ar np is implicit in 
the BA. That is, expanding Eq. pi at q — »• yields 
F n (q)/F p (q) 1 - (r n + r p )(r n -r p )<jr/6, which is driven 
by r n — r p (r n + r p ~ 11.1 fm changes by less than 3% in 
the models). Though Coulomb distortions correct A pv by 
more than 30-40%, the correlation prevails in the DWBA 
result. One sees in Fig. [2] that any nuclear model accu- 
rately calibrated to masses and charge radii nearly falls 
on the best-fit line and that the confidence band of the 
regression is very narrow. Looking at Fig. [I] it can be 
realized that different models, similarly successful for the 
well-known observables, can give the same A pv with dif- 
ferent r n (cf. MSkA, BCP, and SkM*; Sk-Rs, Ska, and 
FSUGold; SkI5 and G2), but almost the same L\r np are 
obtained with these forces. That the prediction band 
of the regression is wider horizontally in Fig. [I] than in 
Fig. [2] points to the same fact. Thus, one expects more 
accurate estimates of neutron observables using the cor- 
relation of Fig. [2] Having found Ar np , one can get r n 
by unfolding the finite size of the proton charge from the 
accurate 208 Pb charge radius [12]. We note that our anal- 
ysis allows one to deduce Ar np and r n from A pv without 
assuming any particular shape for the nucleon density 
profiles. Altogether, we believe our results firmly back 
the commissioning of an improved PREX run where A pv 
can be measured more accurately. The present method 
will permit to retain in Ar np and r n most of the ex- 
periment's accuracy. As recently proposed [26], if r n is 
first precisely known, then a second measurement can be 
made at higher energy to constrain the surface thickness 
of the neutron density of 208 Pb. 

The correlation of A pv with Ar np is universal in the 
realm of mean field theory as it is based on widely dif- 
ferent nuclear functionals. It is of interest to get further 



indications on it by looking at existing experiments. The 
208 pk neu t r0 n densities found via proton elastic scatter- 
ing at 0.8 GeV in [2] and 0.3 GeV in [3] were both deduced 
from the data in a way consistent with the experimen- 
tal charge density of 208 Pb (known by electron elastic 
scattering). We computed A pv using the neutron and 
charge densities quoted in these works and plotted the 
results in Fig. [2] against the central Ar np value of each 
experiment (0.14 fm in [2] and 0.21 fm in [3 ). We did 
the same with the data deduced from the antiprotonic 
208 Pb atom [5] (now using the Fermi nucleon densities 
of Table VI of [5 ). It is seen that the theoretical cor- 
relation of the models nicely agrees with these points. 
Our test value A pv = 0.715 ppm of 3% accuracy from 
PREX would give Ar np as 0.195 ± 0.057 fm (see Fig. [2). 
As reviewed in [IT], we may recall that the recent con- 
straints from strong probes, isospin diffusion, and pygmy 
dipole resonances favor a range 0.15-0.22 fm for the cen- 
tral value of Ar np ( 208 Pb). Recent informations on the 
nuclear equation of state derived from observed masses 
and radii of neutron stars suggest a similar range 0.14- 
0.20 fm [21 [28]. 

Finally, we analyze how PREX can constrain the den- 
sity dependence of the nuclear symmetry energy E sym (p) 
around normal density po> which is characterized by the 
slope coefficient L = 3po dE sym (p)/dp\ Po in the literature 
[TTHSTj . A larger L value implies a higher pressure in 
neutron matter and a thicker neutron skin in 208 Pb. In- 
terest in L permeates many areas of active research, such 
as the structure and the reactions of neutron-rich nuclei 
[T5U2T] . the physics of neutron stars [22-24 , and events 
like giant flares [29] and gravitational radiation from neu- 
tron stars [30]. The available empirical estimates span 
a rather loose range 30 < L < 110 MeV, with the re- 
cent constraints seemingly agreeing on a value around 
L ~ 60 MeV with ±25 MeV spread [T7H2T]. A micro- 
scopic calculation with realistic nucleon-nucleon poten- 
tials and three-body forces predicts L = 66.5 MeV [31]. 
Figure [3] displays the correlation between Ar np ( 208 Pb) 
and L [171419] in the present analysis. Imposing the previ- 
ous constraint Ar np = 0.195±0.057 fm yields L = 64±39 
MeV. While the central value depends on our test as- 
sumption A pv = 0.715 ppm, the spread following from 
a determination of A pv to 3% accuracy, essentially does 
not. Then, we have to conclude that a 3% accuracy in 
A pv sets modest constraints on L, implying that some 
of the expectations that this measurement will constrain 
L precisely may have to be revised to some extent. To 
narrow down L, though demanding more experimental 
effort, a ~1% measurement of A pv should be sought ul- 
timately in PREX. Our approach can support it to yield 
a new accuracy near SAr np ~ 0.02 fm and 5L ~ 10 MeV 
well below any previous constraint. Moreover, PREX is 
unique in that the central value of Ar np and L follows 
from a probe largely free of strong force uncertainties. 

In summary, PREX ought to be instrumental to pave 
the way for electroweak studies of neutron densities in 
heavy nuclei [9] [lOl [26] . To accurately extract the neu- 
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FIG. 3: Neutron skin of Pb against slope of the symmetry 
energy. The linear fit is Ar np = 0.101 + 0.00147L. A sample 
test constraint from a 3% accuracy in A pv is drawn. 

tron radius and skin of 208 Pb from the experiment re- 
quires a precise connection between the parity- violating 
asymmetry A pv and these properties. We investigated 



parity-violating electron scattering in nuclear models 
constrained by available laboratory data to support this 
extraction without specific assumptions on the shape of 
the nucleon densities. We demonstrated a linear correla- 
tion, universal in the mean field framework, between A pv 
and Ar np that has very small scatter. Because of its high 
quality, it will not spoil the experimental accuracy even 
in improved measurements of A pv . With a 1% measure- 
ment of A pv it can allow to constrain the slope L of the 
symmetry energy to near a novel 10 MeV level. A mostly 
model-independent determination of Ar np of 208 Pb and 
L should have enduring impact on a variety of fields, 
including atomic parity nonconservation and low-energy 
tests of the Standard Model [U [32] . 
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